Gaussian paradox and clustering in intermittent turbulent signals 
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A relation between intermittency and clustering phenomena in velocity field has been revealed for 
homogeneous fluid turbulence. It is described how the intermittency exponent can be split into sum 
of two other exponents. One of these exponents (cluster-exponent) characterizes clustering of the 
'zero'-crossing points in nearly Gaussian velocity field and another exponent is related to the tails 
of the velocity probability distribution. The cluster-exponent is uniquely determined by the energy 
spectrum of the nearly Gaussian velocity field and entire dependence of the intermittency exponent 
on Reynolds number is determined by the cluster-exponent. 

PACS numbers: 47.27.-i, 47.27.Gs 



Relation between intermittency of dissipation field 
01 and clustering in turbulent velocity field 0,0 
is still very obscure. Taking into account that turbulent 
velocity itself is nearly Gaussian ,6j the clustering 
phenomenon in the velocity field should be also a 
'Gaussian' one, while the intermittency phenomenon is 
usually associated just with non-Gaussian properties 
of the velocity derivatives 0,0] ,0- On the other 
hand, it is expected that high frequency events in the 
velocity field should provide significant contribution 
to the turbulent dissipation and, especially, to its 
high order moments 0,0,0- Therefore, to find 
relationship between the two phenomena is rather a 
non-trivial task that demands special tools. In our 
recent papers 0,0 we introduced a cluster-exponent to 
describe quantitatively clustering in turbulent velocity 
field and we found an empirical relationship between 
the cluster- and intermittency exponents. In present 
paper we will describe how the intermittency exponent 
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FIG. I: Standard deviation of the running density fluctua- 
tions against r for for the velocity fluctuations measured in 
the wind-tunnel [Tl| (Fig. fa) and for their Gaussian spectral 
simulation (Fig. lb). 
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FIG. 2: As in Fig . 1 but for atmospheric surface layer at 
Re\ — 19000 [l2j| (Fig. 2a) and for its Gaussian spectral 
simulation (Fig. 2b). 



can be split into sum of two other exponents (Eq. 6). 
One of these exponents is the cluster-exponent and 
another exponent is related to the tails of the velocity 
probability distribution. We also will show that the 
cluster-exponent is uniquely determined by the energy 
spectrum of the nearly Gaussian velocity field. This 
simple splitting of the intermittency exponent on the 
Gaussian and 'tail' components can be considered as a 
solution of the problem which was formulated above. 
Moreover, for finite Reynolds numbers Re\ the energy 
spectrum depends on Re\ and, therefore, the cluster- 
exponents depends on Re\, while the 'tail' component 
of the intermittency exponent is independent on Re\. 
Therefore, the entire dependence of the intermittency 
exponent on Re\ is uniquely determined by its Gaussian 
component (through the clustering phenomenon) that 
we will call Gaussian paradox. 

Let us count the number of 'zero'-crossing points of 
the signal in a time interval r and consider their running 
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density n T . Let us denote fluctuations of the running 
density as 6n T = n T — (n T ), where the brackets mean 
the average over long times. We will be interested in the 
scaling variation of the standard deviation of the running 
density fluctuations (Sn 2 ) 1 ' 2 with r 
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For white noise it can be derived analytically 0,^3] that 
a n = 1/2 (see also 0). In Figure la and 2a we show 
calculations of the standard deviation for a turbulent ve- 
locity signal obtained in the wind-tunnel experiment [Tl| 
and for a velocity signal obtained in an atmospheric ex- 
periment ^3 respectively. The straight lines are drawn 
in the figures to indicate scaling (1). One can see two 
scaling intervals in Fig. 1. The left scaling interval cov- 
ers both dissipative and inertial ranges, while the right 
scaling interval covers scales larger then the integral scale 
of the flow. While the right scaling interval is rather triv- 
ial (with a n = 1/2, i.e. without clustering), the scaling 
in the left interval (with a n < 1/2) indicates clustering 
of the high frequency fluctuations. The cluster-exponent 
a n decreases with increase of Re\, that means increas- 
ing of the clustering (as it was expected from qualitative 
observations). 

To describe intermittency a running average of fluctu- 
ations of the dissipation rate e(t) of the velocity signal 
u{t) can be used [L| (see also a discussion in [l4() 
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where 



(2) 



(3), 



and (Se 2 .) 1 / 2 is standard deviation of e T . Figure 3 shows 
an example of the scaling (2) providing the intermittency 
exponent a E for the wind-tunnel data (see |l3j| for more 
details and Fig. 5 for the statistical errors). 

Events with high concentration of the zero-crossing 
points in the intermittent turbulent velocity signal pro- 
vide main contribution to the (Se 2 ) due to high concen- 
tration of the statistically significant local maximums in 
these events ||- Therefore, one can make use of sta- 
tistical version of the theorem 'about mean' in order to 
estimate the standard deviation 



2x1/2 



(4) 



where {max : u 2 } T is maximum of u(t) 2 in the interval of 
length r. Value of (<5{max : u 2 } 2 ) 1 ^ 2 should grow with 
Re\ (due to increasing of (u 2 )) and, statistically, with 
the length of interval r (the later growth is expected to 
be self-similar), i.e. 



(6{max : u 2 } 2 T ) 1/2 = C(Re x ) • r Q " 



(5) 



where constant C(Re\) is a monotonically increasing 
function of Re\, and the statistical ('tail') exponent 
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FIG. 3: Intermittency exponent as slope of the straight line 
for the longitudinal velocity signal obtained in the wind tunnel 
experiment [lj] at Re\ = 1029. 



10 
9 
8 
7 

6 

5 - 
4 - 



a. 2 



wind tunnel 
Re l =1029 




10 1 



10' 



10' 



10 4 



FIG. 4: Statistical ('tail') exponent a m ~ 0.1 as slope of the 
straight line for the longitudinal velocity signal obtained in 
the wind tunnel experiment [Tlj at Re\ = 1029. 



a m > is independent on Re\. Fig. 4 shows an example 
of scaling (5) providing the statistical ('tail') exponent 
a m ~ 0.10 ±0.01. Substituting (5) into (4) and taking 
into account (1) and (2) we can infer a relation between 
the scaling exponents 

a e = a n - a m (6) 

Unlike the exponent a™ the exponent a n depends on Re\ 
(see Fig. 2 and Ref. (4J). Therefore, we can learn from 
(6) that just the clustering of the high frequency velocity 
fluctuations (a„ in (6)) is responsible for dependence of 
the intermittency exponent a £ on the Re\. 

It is naturally consider the exponent functions through 
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FIG. 5: The cluster-exponent: a n (circles), against 1/ln.ReA 
for velocity signal at different values of Reynolds number 
(200 < Rex < 20000). Scaling exponent a e (triangles) for 
fluctuations of the dissipation rate of the velocity signals. 

lnR\: a = a (lni?e A ) B,[ll,[ll (see also Appendix). 
Following to the general idea of Ref. 0] (see also 3) let 
us expend this function for large Rex in a power series: 



a(lnRe\) — ev(oo) + 
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(7) 



In Figure 5 we show the values of a n (circles) and 
a £ (triangles) calculated for the velocity signals against 
1/ In .Re a for 200 < Re x < 20000 (the data are taken from 
the wind-tunnel experiment from the atmospheric 
experiment |12| , from a mixing layer and an atmospheric 
surface experiments 01). The straight lines (the best 
fit) indicate approximation with the two first terms of 
the power series expansion (7) 
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This means, in particular, that 
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The closeness of the constants a\ ~ 3/2 (the slopes of the 
straight lines in Fig. 5) in the approximations (8) of the 



cluster-exponent a n and of the intermittency exponent 
a £ confirms the relationship (6) with a rn ~ 0.1 (cf Figs. 
4 and 5). 

Thus one can see that entire dependence of the in- 
termittency exponent a £ on Rex is uniquely determined 
by dependence of the cluster-exponent a n on Rex- On 
the other hand, the cluster-exponent a n itself is uniquely 
determined by energy spectrum of the velocity signal at 
suggestion that the velocity field is nearly Gaussian. This 
means that energy spectrum uniquely determines depen- 
dence of the intermittency exponent on the Reynolds 
number Rex- From conventional point of view on in- 
termittency it seems as a paradox. 

We use spectral simulation to illustrate the uniquely 
relationship between energy spectrum and cluster- 
exponent for the Gaussian signals; i.e. we generate Gaus- 
sian stochastic signal with energy spectrum given as data 
(the spectral data is taken from the original velocity sig- 
nals). Obviously, the pure Gaussian signal obtained in 
such simulation does not possesses the turbulent inter- 
mittency properties resulting in the estimate (4). Figures 
lb and 2b show, as examples, results of such simulation 
for the wind-tunnel velocity data (Fig. lb) and for the 
atmospheric surface layer (Fig. 2b). Namely, Figures la 
and 2a show results of the calculations for the original 
velocity signals and Figs, lb, and 2b show corresponding 
results obtained for the spectral simulations of the signals 
(to be sure, for the both Gaussian simulations a £ ~ 0.5 
unlike of ct n ). 

It is interesting to note that for isotropic turbulence 
with prominent incrtial (Kolmogorov) range the only 
small-scale edge of the inertial range (this scale depends 
on Rex) determines value of the cluster-exponent a n 
corresponding to the spectrum. I.e. role of the viscous 
dissipation in this respect is reduced to the effective 
small-scale 'cut-off' effect and just dependence of the 
'cut-off' scale on Rex determines dependence of the 
intermittency exponent a £ on Rex through equation (6) 
for such flows. 
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FIG. 6: The vortex ring and its velocity. 
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(see figure 6). 

If, for instance, a straight line-vortex develops a kink 
with a radius of curvature A/2, then self-induction gen- 
erates a velocity perpendicular to the plane of the kink. 
This velocity can be also calculated using (A.l). One 
can guess that in a turbulent environment, the most un- 
stable mode of a vortex tube with a thin core of length 
L (integral scale) and radius r\ (Kolmogorov or viscous 
scale), will be of the order A: Taylor-microscale |(|.[l8j. 
Then, the characteristic scale of velocity of the mode with 
the space scale A can be estimated with help of equation 
(A.l). Noting that the Taylor-microscale Reynolds num- 
ber is defined as Q 
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(A.2) 



where vq is the root-mean-square value of a component 
of velocity. It is clear that the velocity that is more rel- 
evant (at least for the processes related to the vortex 
instabilities) for the space scale A is not vq but v given 
by (A.l). Therefore, corresponding effective Reynolds 
number should be obtained by the renormalization of the 
characteristic velocity in (A.2), as 



A v 



2nv \2n 



(A3) 



It can be readily shown from the definition that 
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I. APPENDIX 

Thin vortex tubes (or filaments) are the most promi- 
nent hydrodynamical elements of turbulent flows |17| . 
The filaments are unstable in 3-dimcnsional space. In our 
recent paper Q we investigated this instability in order 
to renormalize the Taylor microscale Reynolds number. 

A straight line- vortex can readily develop a kink prop- 
agating along the filament with a constant speed. To 
estimate the velocity of propagation of such a kink let us 
first recall the properties of a ring vortex [T^]. Its speed 
v is related to its diameter A and strength V through 

v = K-c ] n-(£-)> ( A1 ) 
2?rA \2n) 

where n is the radius of the core of the ring and A/277 ^> 1 



where const — 15 1 / 4 ~ 2. Hence 

Re e x ff ~ -^Hte\) (4-5) 

47TZ/ 

The strength T can be estimated as 

T - 2iru v r) (A6) 

where v v — u/rj is the velocity scale for the Kolmogorov 
(or viscous) space scale n 0. Substituting (A. 6) into 
(A. 5) we obtain 

Re e / f ~ ln(i?e A ). (A7) 

Thus, for turbulence processes determined by the vortex 
instabilities the relevant dimensionless characteristic is 
hx(Re\) rather than Re\ (cf Eq. (7) and Fig. 5). 



